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Abstract 

Nuclear reactions govern major aspects of the chemical evolution of galaxies and stars. Analytic study 
of the reaction rates and reaction probability integrals is attempted here. Exact expressions for the reaction 
rates and reaction probability integrals for nuclear reactions in the cases of nonresonant, modified nonres- 
onant, screened nonresonant and resonant cases are given. These are expressed in terms of H-functions, 
G-functions and in computable series forms. Computational aspects are also discussed. 

1. Introduction 

Nuclear reactions govern major aspects of the chemical evolution of the universe or, at least, its building 
blocks: galaxies and stars. A proper understanding of the nuclear reactions that are going on in hot cosmic 
plasmas, and those in the laboratories as well, requires a sound theory of nuclear-reaction dynamics (Brown 
and Jarmie 1990). The nuclear reaction rate is the central quantity connecting the theoretical models of 
universal, galactic and stellar evolution and the nucleosynthesizing nuclear reactions, which can be studied 
to some extent in nuclear-physics laboratories. Compilations of reaction rates and uncertainties (analytic 
expressions and tabulated temperature-dependent values) and astrophysical S-factors (analytic expressions 
and tabulat ed energy-dependent values) for ch arged-particle induced reactions are available on the World 



Wide Web (http://csa5.lbl.gov/chu/astro.html). The rate of reacting particles i and j, in the case of 



nonrclativistic nuclear reactions taking place in a nondegenerate environment, is usually expressed as 



(1.1) 



where rii and nj denote the particle number densities of the reacting particles i and j, fi = „™^^ is the 
reduced mass of the reacting particles, T is the temperature, k is the Boltzmann constant, a(E) is the 
reaction cross section and v = (2E/fi)i is the relative velocity. Thus < av > is the reaction probability 
integral, that is, the probability per unit time that two particles, confined to a unit volume, will react with 
each other. The basic assumptions made in (1.1) are that the reacting particles i and j have isotropic 
Maxwell-Boltzmann kinetic-energy distributions and that the kinematic of the reaction can be treated in the 
center-of-mass system, see Fowler (1984), Mathai and Haubold (1988), and Imshennik (1990). 

For nonresonant nuclear reactions between nuclei of charges Zi and Zj at low energies (below the Coulomb 
barrier), the reaction cross section has the form 
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with 
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where f](E) is the Sommerfeld parameter, h is Planck's quantum of action and e is the quantum of electric 
charge. For a slowly varying cross-section factor S(E) we may expand 

S{E)-S{0) + —E+--^ r E . 

Thus 
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Then the integral to be evaluated is 
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"„(*) = / y'e-ve-'v 5 dy. (1.2) 
Jo 

We will consider a general integral of the following form: 
Nonresonant case: 

h= y u e- ay e- zy ~"dy, a > 0, z > 0, p > 0. (1.3) 



It may be a stringent assumption to consider a thermonuclear fusion plasma as being in thermodynamic 
equilibrium. If there is a cut off of the high energy tail of the Maxwell-Boltzmann distribution then (1.2) 
becomes the following: 

Nonresonant case with high energy cut-off: 

r d 

I* = / y v e- ay e~ zy , d<oo, a > 0, z > 0, p > 0. (1.4) 
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We consider also ad hoc modifications of the Maxwell-Boltzmann distribution for the evaluation of the 
nonresonant thermonuclear reaction rate, which acts as a depletion of the tail of the distribution function 
(see Kaniadakis et al. 1996) In this case the general integral to be evaluated is the following: 

Nonresonant case with depleted tail: 

h= y»e- ay e- hy e~ zy "dy, 5 > 0, z > 0, p > 0, a > 0. b > 0. (1.5) 
Jo 

Taking into account the electron screening effects for the reacting particles Haubold and Mathai (1986) 
consider the case of the screened nonresonant nuclear reaction rate (see Lapenta and Quarati 1993). In this 
case 



Tin — TliTli 



^2 1 S^(0) 



x / y » e - y e-< y+ &) 'dy. (1.6) 
Jo 

In this case we consider the following general integral: 
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Screened nonresonant case: 



h= y v e- av e- z{ - v+ V~ P dy, t > 0, p > 0, z > 0, a > 0. (1.7) 
Jo 

When the nuclear cross section a(E) in (1.1) has a broad single resonance it can be expressed via the 
parametrized Breit-Wigner formula and then we have, see also Haubold and Mathai (1986a), 
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Thus the general integral to be evaluated in this case is of the following form: 
Resonant case: 

r°° n.Vp-ay-zy p 

^ = / I — xT^-2 A y- (1-9) 

In the resonant case also we can consider a modification of the Maxwell-Boltzmann distribution which 
results in a depletion of the tail. Then the general integral to be evaluated is the following: 

Resonant case with the depleted tail: 

J 6 = / — ( \2~T 2 dj/. (1.10) 
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The aim of the present article is to give new exact analytic representations of the integrals in (1.2) to 
(1.10). Additional representations are available from Mathai and Haubold (1988). 

2. Reduction formulae for the reaction probability integrals 

Writing 



we have 
Now consider 
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where (a) n denotes the Pochhammcr symbol, 

(a) n = a(a + l)...(a + n - 1), (a) = 1, s^O. 



Note that 
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For simplifying I 5 and I 6 we will use the identity 
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and rewrite the single integral as a double integral. That is 
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Now we expand 
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The integral over 2 gives 



Substituting back we have 
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Thus it is seen that all the integrals I± to I 6 can be reduced to the integral I^"\v,a,z, p) for two different 
situations of non-negative as well as negative v. We will evaluate in the next section by using a statistical 
technique. 



4 



3. Evaluation of the integral 

In general, integrals I\ to Iq are quite difficult to evaluate analytically. Here we will use a statistical 
technique. We will evaluate the density of a product of two independently distributed real scalar random 
variables by using two different methods, one procedure leading to the integral that we want to evaluate 
and the other procedure leading to a representation in terms of a known function. Then appealing to the 
uniqueness of the density we have the integral evaluated in terms of a known special function. Let x and y 
be real scalar random variables having the densities 
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where c\ and c 2 are normalizing factors such that 
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Since the variables are assumed to be independently distributed the joint density of x and y is the product 
of f\{x) and f 2 (y)- Let us transform x and y to u = xy and v = x. Then the joint density of u and v, 
denoted by g(u,v), and the marginal density of u, denoted by gi(u), are given by 

g(u,v) = c 1 c 2 u ,/ v- / - 1 e- av e- cv ~'' , c = zu» 

and 
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Hence we have 
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Let us look at the (s — l)-th moment of u. Due to independence 

£(u' -1 ) = [Eix'-^Eiv'- 1 )]. 

But 



Eixs- 1 )^^ f 
Jo 



d 

x°~'e~ 



= ci > — \ — , for d < oo 

' m! s + m 

m—0 

= cia~ s T(s), SR(s) > 0, for d = oo 
where JJ(-) denotes the real par of (•), and 



E{y s - r ) = c 2 

Jo 



= ^-^r(^), for^ + 5)>0. 
P VP/ 

Taking the inverse Mellin transform of E{u s ~ ) we have 
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where L is a suitable contour and i = This contour integral can be written as an H-function, see for 

example Mathai and Saxena (1978). That is, 
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Substituting (3.3) in (3.2) and then comparing with (3.1) we have 
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(3.5) 
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where 9?(z/) > 0,SR(o) > 0,5ft(z) > 0,5ft(p) > 0,0 < u < oo. Note that (3.4) , (3.5) and (3.6) give an 
I^' \v^a 1 z, p) with v negative. If the integral for a positive v is required then we proceed as follows: Take 
fi(x) — c- i x l/ Q~ ax and ji{y) — c±e~ zx " , where c 3 and C4 are the new normalizing constants, and then proceed 
as before. Then we end up with the following representations. 
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where > 0, SR(a) > 0, 3?(z) > 0, 5ft(p) > 0. Thus (3.7), (3.8) and (3.9) cover the case of a positive v 
in When p is real and rational then the H-functions appearing in (3.3) to (3.9) can be reduced to 

G-functions by using the multiplication formulc for gamma functions, namely, 
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For the theory and applications of G-functions see for example Mathai (1993). For p = ^, m,n = 1,2, ... 
reduction to the G-function is available from Mathai and Haubold (1988). 

The parameters of interest in nuclear astrophysics are I^"\v,a,z,p) for a = 1, z > 0, p = |. In this 
case computable representations will be discussed in the next section. 

4. Computable series representations of the reaction rate integrals 

Let us start with (3.9) for p = 5. Then the H-function to be evaluated is the following: 
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by expanding T(2s) with the help of (3.10). Note that (4.1) is a G-function of the type Gop(')i see f° r 
example Mathai (1993). Observe that for v ^ ± ^, A = 0, 1, ... all the poles of the integrand are simple. 
Then (4.1) generates a series corresponding to each gamma in the integrand. Corresponding to T(s) the 
poles are at s — —n, n = 0, 1, ... and the corresponding residue is 
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Hence the sum of the residues gives 
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where o-fM - ) is a hypergeometric series which is convergent for all a and b. Thus (3.9) is a linear function of 
3 such series of o-FVs for v ^ ±-|, n = 0, 1, .... If v is an integer or half-integer then we have one set of poles 
of order 2 each and one set of order one each. Techniques for handling the integral when the integrand has 
higher order poles are given in Mathai (1993). Series representations of (4.1) for all cases of v are given in 
Mathai and Haubold (1988). 



Now let us examine (3.8) for p = ^. In this case the H-function to be evaluated is given by 
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At s = —v — m there is a pole of order one if v ^ A = 0, 1, otherwise it will be a pole of order 2 at 
this point. When it is a pole of order one the residue is given by 
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which is evidently convergent. 
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At s — —n, n = 0, 1, ... the integrand in (4.3) has poles of order one when v ^ A = 0, 1, ... The 



sum of the residues here is given by 
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The term corresponding to this in (3.8) is the following double series: 

d ,yy t^T i i (* 2 Afl n 

By Horn's theorem on convergence, see for example Srivastava and Karlsson (1985, pp. 56-57) this double 
series is evidently convergent for all a, z and d. The third term of (3.8) as well as all terms for other cases of 
v can be easily seen to give convergent series for all < d < oo, a > 0, z > 0. Similar arguments hold good 
for (3.5) and (3.6) also. 



Let us examine I 3 . We have expressed I 3 in terms of in (2.1). That is, 
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For p — i a typical term in this behaves like a 0-F2 given in (4.2). A typical term will be a constant 
multiple of a double series of the form 
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This by Horn's theorem is convergent for all (x, y)\0 < x < 00, < y < 00 where x — b and y — Hence 
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the series representation of -Z3 is convergent for all a, b, z and <Hor^^±^, A = 0, 1, .... When v is an integer 
or half-integer then the series representation will contain psi functions and logarithmic terms but from the 
structure in (4.4) one can see that the scries will be convergent. 

Now let us examine the complicated form coming from I§ or from (2.4). Here we have the factor 
I^°\v + mi, a, z, p) with toi > 0. For the case p — ^ we have from (3.9) and (4.1) 
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Writing (4.5) as a G-function, see Mathai (1993), we have 



1 



r (s) r s+o r(i/ + mi + l + s) — ds (4.5) 



r (oo) 



v + mi,a,b,- 



" b 0,3 







az 2 




4 


0,l,i/+roi+l_ 



(4.6) 



The behavior of this G-function for small and large values of ^J- is available from Mathai and Saxena (1973, 
p. 307) or Luke (1969, pp. 178-180). For small values of this G-function behaves like unity and hence 



behaves like a -j-f— ■ For large values of the I^" 1 in (4.6) behaves 
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oo 2 m 

EE 

m— mi— 



2m 

mi 



2\ m 



1 cr 



r 21 m / 2 \ 7 







-1 










for 












5 





< 1 



or for - — \g\ < c < - + \g\ which is equivalent to < c < 1 + \g\ for \g\ > 1, when a = 1, and 
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< 1 or for 



2a[c- < z < 2a[c+ |.g|]*. 

For example, combined with the condition for small values, we have < z < 2(2 + \g\) for a = 1, \g\ > 2. 

For other values of p also the procedure remains the same. We may observe that instead of the I5 in 
(1.9) we can also consider a more general integral of the type 
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and a corresponding expression for (2.4) is obtained. Convergence conditions can be checked exactly the 
same way as in the case of d — 1 = 0. 

Iq can be reduced to a form corresponding to (2.4). Convergence conditions can be checked by converting 
the series form into one of the standard triple series discussed in Srivastava and Karlsson (1985) or by using 
the general procedure for a multiple series or by writing the kernel function as a G-function, as described 
above, and then checking the behavior for large and small values of the argument of this G-function. Note 
that I7 can be expressed in terms of as follows: 
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Writing in terms of a G-function and then checking the behavior of the G-function for small and large 
values of the argument one can verify the existence of I7. For small values, the G-function in I 2 behaves like 
unity and then for p — \ 



,2m-mi 
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for [l — < 1 or i — \g\ < c < i + For large values the behavior of the G-function is available 
from (5.2) later on and in this case the conditions for the existence are given in (5.8) later on. 

When p = — , m,n = 1, 2, ... it is easy to see that the H-function in all the integrals I\ to Iq reduce to 
G-functions, of course with more parameters. General series representations of all forms of G-functions are 
available from Mathai (1993). 

5. Computational aspects 

Anderson, Haubold and Mathai (1994) looked into the computational aspects of the integrals I\ to I4. 
Exact computations and graphs are given there for p = |, a = 1 in the cases of I\ for v = 1, 2; I2 for 
v = 0, d = 1, 5; I3 for (v,S,b) = (0,2,0.001), (1,5,1); h for (v,t) = (0,1), (0,5). For large values of z 
one can use the asymptotic forms of the integrals for computational purposes. These can be worked out by 
using the asymptotic form of the G-function. From (3.9) and (4.1) 



r(°°) 



Lr 0,3 



2 




4 


f+l,o,i. 



(5.1) 



But for large values of ^f- we have 



3,0 
0,3 



have 




2 








4 


f+l,0,i. 



^«-<^)e-»m* 



2^ + 1 

2\ —g— 



(5.2) 



By substituting this expression for the G-function in I\ to J7 we get the following forms for large values of 
^ with a = 1, p= |: 
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For a broad overview of the aspects of numerical evaluation of special functions, including available software 



packages for this purpose, see on the World Wide Web: http://math. nist.gov/nest/ 
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